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We examine dark energy models in which a quintessence or a phantom field, <f>, rolls near the vicin- 
ity of a local minimum or maximum, respectively, of its potential V (<j>)- Under the approximation 
that (l/V)(dV/d(j>) <C 1, [although (l/V)(d 2 V/d(f) 2 ) can be large], we derive a general expression for 
the equation of state parameter iu as a function of the scale factor for these models. The dynamics of 
the field depends on the value of (1/V) (d 2 V/d<p 2 ) near the extremum, which describes the potential 
curvature. For quintessence models, when (1/V) (d 2 V / d(f> 2 ^ < 3/4 at the potential minimum, the 
equation of state parameter w(a) evolves monotonically, while for (1/V) (d 2 V/d(j) 2 } > 3/4, w(a) has 
oscillatory behavior. For phantom fields, the dividing line between these two types of behavior is at 
(1/V) (d 2 V/d(f> 2 ) = -3/4. Our analytical expressions agree within 1% with the exact (numerically- 
derived) behavior, for all of the particular cases examined, for both quintessence and phantom fields. 
We present observational constraints on these models. 



I. INTRODUCTION 

It has been known for over a decade [J, [|J that at least 
70% of the energy density in the universe is in the form 
of an exotic, negative-pressure component, called dark 
energy (see Ref. Q for a recent review) . The equation of 
state of this dark component, defined as the ratio of its 
pressure to its density: 



w = pue/p 



DE; 



(1) 



is observationally constrained to be close to —1. In par- 
ticular, constraints on constant w lead typically to w = 
— 1 with ±10% accuracy 0, However, a variety of 

cosmological paradigms, based on scalar fields, attribute 
a dynamical nature to dark energy, leading to a time- 
varying w. The class of models in which the scalar field 
is canonical is dubbed quintessence @, H, H, H3) EH E3 
and has been extensively studied in the literature. An al- 
ternative approach is phantom dark energy, a component 
for which w < — 1 , which can be realized by using a scalar 
field with a negative kinetic term in its Lagrangian, as 
first proposed by Caldwe ll Il3i . Such models have well- 
known problems [3, [TH, Il6l. Il7j , but nevertheless have 
also been widely studied as potential dark e nergy c andi- 

dates 0, M M M M M M HE M MM, HH- 

Finally, since both quintessence and phantom models 
cannot lead to a dark energy equation-of-state parame- 
ter that crosses — 1 during cosmological evolution (which 
might be the possibility according to observations), the 
simultaneous consideration of both models in a joint sce- 
nario, named quintom, has recently gained significant at- 
tention m m, m m, m, m . 
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Given the considerable freedom that exists in choosing 
the potential function V (<fi) of the scalar field, it is use- 
ful to develop general expressions for the evolution of w 
which cover a wide range of models. A general result of 
this type was derived in (3(| for a class of quintessence 
models in which w is assumed to be always close to — 1 
and the potential is nearly flat. The flatness of the po- 
tential is characterized by two "slow roll conditions" : 



and 



1 dV 
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(3) 
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For these models, it was shown in [36( that the evolution 
of w is described by a unique expression involving only 
the present values of and w. In [3?J this result was 
extended to phantom models satisfying (J3j , and the 
corresponding w dependence was shown to be described 
by the same expression with an overall sign change. Fi- 
nally, the extension to quintom scenario under ©-(El) 
was performed in [38j . where a universal expression for 
w was also extracted, allowing for crossing of w = — 1. 

While the slow-roll conditions ([2|) and ([3]) are sufficient 
to give w ~ — 1 today, they are not necessary. Classes 
of models characterized by the validity of © alone, i.e., 
without ([3]), were considered in [3^ |. corresponding in 
particular to a quintessence field in the vicinity of a local 
maximum of its potential. In this case, there is an extra 
degree of freedom, namely the value of (l/V)(d 2 V/d4> 2 ), 
which describes the curvature of the potential in the 
vicinity of the maximum. As a result, instead of a single 
solution for the evolution of w one obtains a family of so- 
lutions that depend on the present-day values of and 
w and the value of (l/V)(d 2 V/d(f> 2 ) at the maximum of 
the potential. As expected, this family of solutions in- 
cludes the slow-roll solution of [36| as a special case in 



2 



the limit where (l/V){d 2 V/d<f> 2 ) -> 0. In 0, a simi- 
lar result was derived for a phantom field evolving near 
a minimum of its potential. It was found that a unique 
family of solutions, very similar to the one derived in [39j, 
can be used to approximate the behavior of w in these 
models. 

More recently, Chiba [4l| showed that while conditions 
([2]) and ([3]) are sufficient to allow the quintessence ver- 
sion of the slow-roll conditions to be applied, they are 
not necessary. Chiba extended the methodology of Ref. 
[39[ | to provide a more general set of conditions on the 
potential, dropping the assumption that the field is close 
to a local maximum in the potential, while still assuming 
that w is close to —1 throughout the evolution. Inter- 
estingly, the expression derived under these more general 
conditions coincides exactly with the expression derived 
in [3^|, indicating the generality of this result. 

In this work we extend the techniques developed in [39[ 
to the "opposite" of the case considered there. In Ref. 
[39[ | . the quintessence field was considered to roll near a 
local maximum in the potential, so that V" > 0. In this 
paper, we consider the case for which V" < 0. This cor- 
responds to a canonical scalar field rolling close to a min- 
imum of its potential. For completeness, we also examine 
the corresponding phantom model, for which the phan- 
tom field rolls near a maximum of its potential. That 
is, while in [39J and [40J the field rolls away from an un- 
stable potential-extremum, in the present work the field 
rolls towards a stable extremum. (Recall that phantom 
fields roll away from local minima and are attracted by 
local maxima). As we will see, this investigation yields 
a more complicated dynamics for the scalar field, which 
can include oscillations in addition to slow roll behavior. 
For simplicity, we assume that condition ^ is satisfied, 
and that the field evolves near an extremum in the poten- 
tial, rather than making the more general assumptions of 
Ref. [4l|. The result can be generalized, in a straighfor- 
ward way, following the development of Ref. [4l|. We 
note further than our main result is mentioned in Ref. 
[4lj . although it is not investigated in detail there. 

The plan of our paper is as follows: In section HT1 wc 
construct the model of scalar field evolution near a sta- 
ble potential extremum, and we derive the expression for 
w(a). In section IIIII we apply our formula to various 
potentials, comparing the results to the exact evolution 
arising by numerical integration. We then constrain this 
general family of models with SNIa observations. Our 
results are summarized in section HVl 



II. FIELD EVOLUTION NEAR A STABLE 
POTENTIAL-EXTREMUM 

We are interested in models where the scalar field 
evolves near a stable potential-extremum, and thus we 
consider a minimally coupled scalar field 4> in a potential 
V (</)), where the field <j> can be either a canonical or a 
phantom one. In the following we introduce the usual e- 



parameter, acquiring the value +1 for the canonical case, 
and —1 for the phantom one. Doing so we can present 
our expressions in a general way. 

The Euler-Lagrange equation of motion of the field 
reads: 



■ + 3H(j> + e— = 0, 



(4) 



where a is the scale factor and H = a/a is the expansion 
rate. Dots denote derivatives with respect to time and 
primes denote derivatives with respect to the field 4>. In 
a flat universe, the expansion rate is linked to the total 
density px via the first Friedman equation (in units where 
8ttG = 1) as 



H 2 = p T /3. 



(5) 



Additionally, the evolution of the scale factor is given by 
a 1 



G 



(PT + Pi) ■ 



where pt is the total pressure. 

Following [39l | we perform the transformation 

4>{t) = u{t)/a{tf'\ 
and therefore Eq. (Q| becomes 

it + jp T u + ea 3/2 V (u/a 3/2 ^j = 0. 



(6) 



(7) 



(8) 



We consider a universe consisting of pressureless mat- 
ter and a scalar field, where the scalar field plays the role 
of the dark energy. In order to realistically mimic the 
observed dark energy, the scalar field must have w close 
to —1 and its energy density must be roughly constant. 
The total pressure pt is then simply given bypx ~ ~/°</>o, 
where p^o is the present day density of the dark energy. 
(In what follows, a subscript always indicates a present 
day value). Under this approximation, Eq. 1(5)) becomes: 



0. 



(9) 



Since we are interested in applying Eq. Q to a scalar 
field evolving near a local potential-extremum at </>*, for 
any <f> close we use the expansion: 

V (4>) — V (<£*)+(l/2) V" (0*) (0 - <f>*f+0 ((0 - <£») 3 ) . 

(10) 

Substituting the above relation into Eq. ((9]), and im- 
posing V (0* ) = p^o we obtain the following differential 
equation for the field evolution: 

u - [(3/4) V - eV" (0,)] u = 0. (11) 

This is essentially the same equation derived previously 
in Ref. ||. 

In this work we examine cosmological evolution near 
a stable potential-extremum, that is either a canonical 
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field (e > 0) near a local minimum (V"(4>*) > 0), or a 
phantom field (e < 0) near a local maximum (V" < 
0). Thus, in both cases of interest eV"{(j>*) > and 
therefore the following analysis can be performed in a 
unified way. 
Defining 



k = y/(3/4)V(<l>.)-eV" (0*), (12) 

we obtain the general solution of Eq. (| 1 1 [) as 

u = Asinh(fct) + Bcosh(fci) , (13) 

where A and B are arbitrary constants. Note that this 
solution holds for k being either real or imaginary, where 

I 



in the latter case the hyperbolic functions are straight- 
forwardly replaced by trigonometric ones. Hence, we can 
simply use the earlier results of Ref . [39j . 
Defining 



t A = 2/ v /3^ = 2/y/3V (</>*), 

and taking <j) (t = 0) = <pi , we obtain 

(pi sinh (kt) 
kt\ sinh (t/t\) ' 



(14) 



(15) 



and the equation of state parameter, w(a), is given by 
(see Ref. [39j for the details): 



1 + w(a) = (1 + wo)a 



-3 



^Tl^oktA cosh [kt (a)] — yj (1 — ^<>o)a 3 + ^</>o sinh [fci (a)] 



1 2 



[y^n^ofe^A cosh(fcio) — sinh(fcto)] 
I 



(16) 



where wq is the present-day value of w(a), and t(a) and in favor of u>q which is obviously smaller than —1 for 



to can be derived from: 



t(a) = t\ sinh 1 



1 — f^o 



phantom while it is larger than —1 for the quintessence 
case. 

Following Ref. [3{|, we now introduce the constant 

(17) 

K = kt\. In terms of the potential, K can be written as 



and 



to = t\ tanh 



1 f-\A^o) 



4 V" ( 



(18) 



3V(0* 



(19) 



In terms of K we can finally express the evolution of w 
Note that in (fT!)]) the e-parameter has been simplified in the following form: 

I 



, M h , \ 3 (K-i) {(Fja) + l) K (K - Fjaj) + (F(a) - l) K (K + F(a))} 
+ w(a) = (1 + w )a y ' T~P) 775 775 77? 



r 



(20) 



where F(a) is given by 



F(a) 



i + (n. 



-1 
<t>o 



Da- 3 . 



(21) 



(Note that F(a) = 1/ wTl^(a), where ft,p(a) is the value of 
il,p as a function of redshift, so that F(a = 1) = ^ •) 
Expression (|20[) is identical in form to the expression 
for w(a) for a quintessence field near a local maximum 

38 | or phantom near a local minimum [40] ] (and see Ref. 

4l| for the derivation of this expression for more general 
cases). However, the crucial difference lies in the defini- 
tion of K (equation [19)1 . In the cases considered in Refs. 



[39l | and |4(j, K was always real since eV" ((/>*) < (in 
the present language). In the case we are considering 
here, for which sV"(4>*) > 0, K can be real or imagi- 
nary (see also [4l|), corresponding to eV" /V < 3/4 or 
eV"/V > 3/4, respectively. Further, for the special case 
of K = 0, equations (fTS"]) and (|2"0"|) are no longer valid; 
instead we have 



(22) 



t\ sinh (t/t\) 



and 



4 



1 + w(a) = (1 + wo)a 



-3 



,/(! - f^ )a 3 + O^o sinh 



-1 / fa^ 



(23) 



Now consider the behavior of the scalar field for the 
three cases of interest, K 2 > 0, K = 0, and K 2 < 0. 
When K 2 > 0, the functional form of equation (f2"D|) 
for quintessence evolution near a potential minimum (or 
phantom evolution near a local maximum) is identical 
to the form for the evolution of quintessence near a 
maximum (or phantom near a minimum) given in Refs. 
[39l [40( | . or for the evolution described in Ref. [4l| . 
However, the crucial difference is that in the previously- 
considered cases, we have K > 1, while for the case con- 
sidered here, we have instead < K < 1 (the special 
case K = 1 is discussed in Ref. [HI). Therefore, al- 
though qualitatively the behavior of w(a) is similar to 



that in Refs. [IEEEEl|]> i.e. w(a) decreases monotoni- 
cally, there are significant quantitative differences. 

When K = 0, we can no longer use equation ([2H)l . but 
equation (|23[) gives an evolution for w(a) that is qualita- 
tively similar to the K 2 > case, i.e., a slow monotonic 
decrease in w. 

Finally, when K 2 < 0, w(a) is oscillatory. While equa- 
tion ([20]) is formally valid (and gives a real expression 
for w(a)) in this case, the oscillatory behavior becomes 
more transparent by writing K = it\j, where 7 is real, 
and simplifying Eq. (|f 6p to a more intuitive form (see 
also Ref. [Il[): 



1 + w(a) = (1 + w )a 



-3 



/f^oT^A cos [jt (a)] - >/(l-fi0o)a 3 + Sl0 O sin [7* (a)] 



/ H^7t A cos(7i ) - sin (7*0)]' 



(24) 



The behavior of cf>(a) for several different values of K is 
illustrated in Fig. fT]). The quantity if can be physically 
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FIG. 1: The evolution of the scalar field <f> as a function of 
the scale factor a, for the indicated values of K, as defined in 
Eq. (HI}. 

interpreted as a measure of the sharpness of the curvature 
of the potential at its extremum. From the definition of 
K (Eq. Q35J we note that in order for K to be real, \ V"\/V 
is small, implying a "flat" potential. The field evolves 
in a slow, friction-dominated manner (with the Hubble 



parameter acting as the friction coefficient according to 
Eq. 2]), asymptotically coming to rest at the potential 
extremum (minimum for quintessence and maximum for 
phantom). On the other hand, an imaginary K requires a 
large |V"|/V, implying a sharp curvature of the potential 
at the extremum, which allows for oscillations since the 
friction term can be overcome. 

It is important to note that we fix the minimum value 
of the potential to be at the level of the cosmological 
constant, so the evolution is still potential energy dom- 
inated (i.e., w ~ — f) as the field oscillates around the 
extremum. These oscillatory solutions are therefore very 
different from the ones usually considered in oscillating 
dark energy models [H, EE El, EE EE E3]> where the 
potential minimum is fixed at V = 0, and w oscillates 
between ±1. 



III. COMPARISON TO EXACT SOLUTIONS 

In this section, we compare our analytic expression for 
the evolution of w to the numerically computed exact 
evolution for a few different models. In each case we have 
a perfect-fluid dark matter and a dark-energy component 
attributed to a quintessence (e = +1 ) or a phantom (s = 
— I) field (p. For the quintessence case, we consider three 
different potentials which have local minima, and we use 
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the corresponding inverted potentials for the phantom 
case. Our purpose is not to propose these toy models 
as specific possibilities for the dark energy, but rather to 
test the accuracy of our approximation against a variety 
of different possibilities. 

The PNGB model [H| , has a potential given by 



V(ct>) =p^ + eM A [1- cos (</»//)]. 



(25) 



where M and / are constants. (For recent discussions of 
the PNGB model in the context of dark energy, see, e.g., 
Refs. [49|, HH and references therein). Other mod- 
els with a local potential minimum include the Gaussian 
potential, 



V{fa) = P4>0 + sM 4 
and the quadratic potential 



1 



V(4>) = pto + £V 2 <j> 2 . 



(26) 



(27) 



where a and V2 are constants. We set initial conditions 
deep within the matter-dominated regime. The value of 
the potential at the extremum, p^Q, is chosen to be equal 
to the energy of the cosmological constant. The initial 
velocity of the field is taken to be zero. 

As discussed above, our formalism applies to models 
for which Q is satisfied, but ([3]) is not. The initial value 
of the field fa determines the accuracy of the first slow- 
roll condition ([2]) at the initial time. To push our formal- 
ism to its limits, we choose fa to a value for which 



1 dV 



0[1) 



For smaller fa of course, the agreement is better. 

We have examined three types of curvature of the po- 
tential at the extremum, characterized by K. The "flat" 
regime is represented by < K 2 < 1, where w evolves in 
a slow monotic manner, and K 2 = represents a limiting 
case of this behavior. The "curved" regime is represented 
by — 00 < K 2 < in which w eventually oscillates. In 
this regime, we have considered the cases of K 2 = — 10 
and K 2 — —20 for each model. 

In Figs. (HE]), the evolution of w from J33J (for K 2 = 
0) or (pM|) (for K 2 < 0) is shown in comparison to the 
exact evolution for the three different models, and the 
different K— regimes described above. The agreement, in 
all three cases, between our analytic approximation and 
the exact numerical evolution is excellent, with errors 
5w < 0.01 in all cases. This result indicates that our 
expressions for w(a) for our three cases of interest can 
be considered "generic" expressions that apply to a wide 
range of possible quintessence and phantom models. 

Finally, in Figs. ([51 fTU|) . we use our analytic approx- 
imation (Eqs. [231I24]) to construct a x 2 likelihood plot 
for wq and f^o for the three choices of K 2 using 
the recent Type la Supernovae standard candle data 
(ESSENCE+SNLS+HST from @). We have exploited 
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FIG. 2: The evolution of w for quintessence in a quadratic po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 
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FIG. 3: The evolution of w for quintessence in a PNGB po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 
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FIG. 4: The evolution of w for quintessence in a Gaussian po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 



6 




0.2 0.4 0.6 0.8 1 



a 

FIG. 5: The evolution of w for a phantom in a quadratic po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 
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FIG. 6: The evolution of w for a phantom in a PNGB po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 
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FIG. 7: The evolution of w for a phantom in a Gaussian po- 
tential for three different values of K 2 . The solid blue curves 
indicate the exact evolution and the red dashed curves indi- 
cate the analytic prediction. 




FIG. 8: Likelihood plot from SNIa data for the parameters wo 
and fi^o, for quintessence and phantom models with generic 
behavior described by Eq. (|23|l . with K 2 — 0, where K is 
the function of the curvature of the potential at its extremum 
given in Eq. (|19[l . The yellow (light) region is excluded at 
the 2<j level, and the darker (orange) region is excluded at 
the lcr level. Red (darkest) region is not excluded at either 
confidence level. 
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FIG. 9: Likelihood plot from SNIa data for the parameters wo 
and Q^o, for quintessence and phantom models with generic 
behavior described by Eq. (I24|l . with K 2 — —10, where K is 
the function of the curvature of the potential at its extremum 
given in Eq. (|19[1 . The yellow (light) region is excluded at 
the 2<j level, and the darker (orange) region is excluded at 
the la level. Red (darkest) region is not excluded at either 
confidence level. 



the fact that our expressions for w(a) for quintessence 
and phantom models have the same functional form, al- 
lowing us to produce a likelihood plot that is continuous 
across wq = —1. However, it is important to note that 
in these figures, the dashed line at Wq = — 1 divides two 
distinct models. Clearly, these models are not ruled out 
by current supernova data. 
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FIG. 10: Likelihood plot from SNIa data for the parame- 
ters wo and fl^o, for quintessence and phantom models with 
generic behavior described by Eq. (|24|) . with K 2 = —20, 
where K is the function of the curvature of the potential at 
its extremum given in Eq. (|19p . The yellow (light) region 
is excluded at the 2<r level, and the darker (orange) region is 
excluded at the la level. Red (darkest) region is not excluded 
at either confidence level. 



to a stable local potential extremum (i.e., a minimum for 
the quintessence and a maximum for the phantom case, 
respectively). Such models provide a mechanism to pro- 
duce a value of w that is close to —1. We have tested our 
expression against the (numerically determined) exact 
evolution for three different representative models and in 
each case it replicated the exact evolution studied with an 
accuracy greater than 99%. A comparison between our 
generic approximation and the observational data indi- 
cates that these models are allowed by SNIa data for a 
variety of values of the potential curvature parameter K 
defined in equation (fF9"| . 

Finally, we note that the case considered here, in which 
V" > for quintessence (and V" < for phantom mod- 
els) leads to a much richer set of behaviors than the 
previously-examined case of V" < for quintessence (and 
V" > for phantom models). In the case examined here, 
we see three very different regimes, depending on the sign 
of K-. 



IV. CONCLUSIONS 

Using techniques previously applied in [H, H3, 53 , we 
have derived a general expression for the evolution of w, 
which is valid for a wide class of quintessence and phan- 
tom dark energy models in which the field is rolling close 
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